Astronomy & Astrophysics manuscript no. Lindetal' final 


© ESO 2009 


June 22, 2009 





Signatures of intrinsic Li depletion and Li-Na anti-correlation in the 
metal-poor globular cluster NGC 6397 * 

K. Lmd\ F. Primas\ C. Charbonnel^'^ , F. Grundahl'*, and M. Asplund^ 

6> . 

J.— s ^ European Southern Observatory (ESO), Karl-Schwarzschild-Strasse 2, 857 48 Garching bei Munchen, Germany 

e-mail: klind@eso.org 
^ Geneva Observatory, 51 chemin des Mailettes, 1290 Sauvemy, Switzerland 

^ Laboratoire d'Astrophysique de Toulouse-Tarbes, CNRS UMR 5572, Universite de Toulouse, 14, Av. E.Belin, F-31400 Toulouse, 
^ ' France 

^ , * Department of Physics & Astronomy, Aarhus University, Ny Munkegade, 8000 Aarhus C, Denmark 

^ Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Strasse 1, 857 41 Garching bei Munchen, Germany 

(N 
(N 



00 



Received 2009 May 19; accepted 2009 June 10 

ABSTRACT 



Context. To alleviate the discrepancy between the prediction of the primordial lithium abundance in the universe and the abundances 
observed in Pop II dwarfs and subgiant stars, it has been suggested that the stars observable today have undergone photospheric 
(~ i , depletion of lithium. 

^ |. Aims. To identify the cause of this depletion, it is important to accurately establish the behaviour of lithium abundance with effective 

I ' temperature and evolutionary phase. Stars in globular clusters are ideal objects for such an abundance analysis, because relative stellar 

O ' parameters can be determined precisely. 

Methods. We have conducted a homogeneous analysis of a very large sample of stars in the metal-poor globular cluster NGC 6397, 
^ . covering all evolutionary phases from below the main sequence turn-off to high up on the red giant branch. Non-LTE Li abundances or 

^ ' abundance upper limits were obtained for all stars, and for a sizeable subset of the targets sodium abundances were also obtained. The 

Na abundances were used to distinguish stars formed out of pristine material from stars formed out of material affected by pollution 
, from a previous generation of more massive stars. 

^ ■ Results. The dwarf, turn-off, and early subgiant stars in our sample form a thin abundance plateau, disrupted in the middle of the 

' subgiant branch by the Li dilution caused by the first dredge-up. A second steep abundance drop is seen at the luminosity of the 

' red giant branch bump. The turn-off stars are more Li-poor, by up to 0. 1 dex, than subgiants that have not yet undergone dredge- 

QQ , up. In addition, hotter dwarfs are slightly more Li-poor than cooler dwarfs, which may be a signature of the so-called Li dip in the 

. cluster, commonly seen among Pop I stars. The feature is however weak. A considerably wide spread in Na abundance confirms that 

• ' NGC 6397 has suffered from intracluster pollution in its infancy and a limited number of Na-enhanced and Li-deficient stars strongly 

contribute to forming a significant anti-correlation between the abundances of Na and Li. It is nevertheless seen that Li abundances 
, are unaffected by relatively high degrees of pollution. Lithium abundance trends with effective temperature and stellar luminosity are 

0^ . compared to predictions from stellar structure models including atomic diffusion and ad-hoc turbulence below the convection zone. 

' We confirm previous findings that some turbulence, with strict limits to its efficiency, is necessary for explaining the observations. 



Key words. Stars: abundances - Stars: Population II - Globular clusters: general - Techniques: spectroscopic - Methods: observa- 
tional - Diffusion 

1. Introduction standard Big Bang nucleosynthesis (BBNS) (ICvburt et al.l2008l) . 

^, , . ^ - . . , , , , , , This is significantly higher than what is found for the "Spite 
The detection of Li in the atmospheres of old and metal-poor j^^^^^,, ^^^-^^ ^ ^^-^^ ^g^^) in the Galactic halo, a well- 
Population I stars has opened up an observational link to the ^^^^^^ ^i abundance plateau consisting of metal-poor halo 
primordial Universe. ^Li is indeed one of four isotopes that ^^^^.^^ ^^^j ^^^ -^^^ ^^^^^^ (g^g^ ^^^^.^ 

were synthesised by nuclear reactions shortly after the Big Bang. „, , , . ■ , ^ ,j , , j 

, . . J- , , J ■ r The plateau abundance in the held has been deter- 

With the most recent determination of the baryon density from • . • i j- i 

the 5-year relea se of WMAP data, ^,h^ = 0.02273 ± 00062 ^" % n .'4 f , . , ° ..'^'"r 7"^" 

(iDmiMeykaLlM, an initial abundance of A^(^Li)/A^(H) = ^^^i) ^ 2.0 -2.4 for stel ar metallicities [Fe/H] = 

< 1/1+0 71 in-io — An ■\ nin . n n^JTl • • Af (-3.5) - (-1.0), With a possible tendency towa rd increas- 

5.24*" 't X 10 or A(Li) = 2.72 + O.OqJ is obtained from : , • 1. j -.u • • . n- / r~n 

ing Li abundance with increasing metallicity (e.g. Ryan et alj 

* Based on data collected at European Southern Observatory (ESO), ^001; Melendez & Ramirez 2004; Charboiinel & Primas 200| 

Paranal, Chile, under program IDs 077.A-0018(A) and 281.D-5028(A), Asplund et al. 2006; Bo nifacio et aU2007a| ; |Hosford et al.i2009| ; 

as well as data collected with the Danish 1 .54 m at European Southern ,Aoki et al.,2009). Differences between the various analyses may 

Observatory (ESO), La Silla. arise from the choice of effective temperature scale and correc- 

** Tables 2, 3, and 4 are only available in electronic form at the tions for non-LTE effects. It is nevertheless clear that a discrep- 

CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via ancy by a factor of 2 - 4 with the cosmological prediction is 



|http://cdsweb.u -strasbg.fr/ unavoidable, suggesting that the stars have undergone the corre- 

' A(Li) = log [jj^) + 12 sponding surface depletion of Li. An important fact is that in all 
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the recent observational analysis (i.e., 21st century) of halo field 
stars, no evidence of significant dispersion has been found along 
the plateau, except for a few stars with str ong Li abnormalities 
(e.g. lRvan et anil999HAsplund et al.ll2006l) . 

So on after the discovery of the Spite plateau. iMichaud et al.l 
(1 19841) predicted Li depletion under the combined action of grav- 
itational settling and weak turbulence in the radiative zones 
of Population II stars. Using sophisticated stellar models that 
treat at omic diffusion and ra diative acceleration from first prin- 
ciples, iRichard et alj (l2005l) illustrate how the comparably low 
Li value found in old, metal-poor stars could be naturally ex- 
plained, assuming rather strict limits to the turbulent transport 
extent and efficiency. However, as of today the nature of the un- 
derlying physical mechanism responsible for turbulence has not 
been definitively identified, the main difficulty for the conjec- 
tured processes (mass loss, rotation-induced mixing, ...) being to 
acc ount for the negligible Li dispersion along the plateau. 

ICharbonnel & PrimasI (l2005i) revisited the literature Li data 
for halo field stars with particular focus on the evolutionary sta- 
tus of the sample stars. This study discovered for the first time 
that the mean Li value appear to be higher for the turn-off (TO) 
and SGB stars than for the dwarfs. This result, together with the 
finding that all halo stars with Li abnormalities (i.e., strong de- 
ficiency or high content) lie on, or originate from, the hot side 
of the plateau, lead the authors to suggest that the most massive 
(i.e., post-MS) of the halo stars still observable today have had 
a different Li history t han their less massive dwarf counterparts. 
ICharbonnel & PrimasI suggested that such a behaviour may be 
the signature of a transport process of chemical elements and 
angular momentum whose efficiency changes on the blue edge 
of the plateau. This behaviour corresponds to that of the gener- 
ation and filter ing of internal gravity waves in both Population 
II and I stars ("Talon & Charbonner2003l ICharbonnel & TalonI 
l2008h . ITalon & Charbonnel (2004) describe how internal grav- 
ity waves coupled with rotation-induced mixing are expected to 
lead to higher Li homogeneity among the plateau dwarf stars 
than among the more massive, slightly evolved stars, a scenario 
which may explain the observational findings by Charbonnel & 
Primas. Speaking in favour of this model is its ability to simul- 
taneously explain the internal solar rotation profile and the time 
evolution of the Li abundance at the sur face of solar- and F- 
type stars as seen in Galactic open clusters (ICharbonnel & TaionI 
l2005l) . 

By turning to stars in metal-poor Galactic globular clus- 
ters, the intrinsic stellar processes involved can be further con- 
strained. This approach presents some obvious advantages. In 
particular, the evolutionary status of the observed stars can be 
determined unambiguously and, presumably, all stars in a clus- 
ter were born with the same metallicity (although the surface 
metallicities observed today may vary between stars in different 
evolution ary phases, due to the effects of atomic diffusion (e.g. 
iKorn et a l. 2007; Lind et al. 2008)). 

NGC 6397 is one of the most well-studied metal-poor glob- 
ular clusters and its Li content has been documented down 
to the magnitude of the cluster TO point in several studies 
(Pasquini & Molaro"1996'; 'Castilho et al."2000'; 'Th evenin et all 
[2001; Bonifacio et al. 2002; Kornetal. 2006, 2007)), however, 
with rather poor number statistics. Bonifacio et al. found A(Li) - 
2.34 + 0.06 for the mean abundance of twelve TO stars, whereas 
Korn et al. found 2.24 + 0.05 for five TO stars and 2.36 ±0.05 for 
two SGB stars. By comparing this observed abundance differ- 
ence of Li (as well as Fe, Ca, Mg, and Ti) t o pred ictions from the 
stellar-structure models by Richard et al.l (l2005l) and references 
therein, Korn et al. empirically constrain the efficiency of slow 



macroscopic motions counteracting atomic diffusion below the 
convective envelope, which should ultimately provide clues on 
the origin of turbulence. Note that some slow macroscopic pro- 
cess is also required to reproduce the observe d morphologies of 
globular cluster colour-magnitude diagrams dVandenBerg et al.l 
120021) . 

Interpreting abundance trends in globular clusters must how- 
ever be done with great caution. Indeed, it is well-known that 
globular-cluster stars present striking anomalies in their light 
element content that ar e not seen am ong their field counter- 
parts (for reviews see e.g. lGratton et al . 2004; Charbonnel 2005|. 
More precisely, C, N, O, Na, Mg, and Al abundances show 
large star-to-star variations within individual clusters. C and 
N, O and Na, and Mg and Al are respectively anti-correlated, 
the abundances of C, O, and Mg being depleted in some stars 
while those of N, Na, and Al are enhanced. Importantly, the 
abundance of Li was found anti-correlated with that of Na and 
coiTelated with that of O in turnoff stars in th e more metal- 
rich globular clusters NGC 6752 and 47 Tuc dPasquini et al.l 
120051 : Bonifacio et al.ll2007bb . These abundance patterns are ex- 
plained by contamination of the star-forming gas by hydrogen- 
processed material ejected by a first generation of short-lived 
massive globular-cluster stars. This increases the Na abundances 
while lowering the O and Li abundances with respect to the 
pristine cluster composition in a second generation long-lived 
low-mass stars. In the framework of the self-enrichment sce- 
nario the Li content of globular-cluster stars is actually an im- 
portant tool to quantify the dilution factor between the ejecta of 
the massive star s responsible for pollution and t he pristine intra- 
cluster matter (Prantzos & Charbonnel 2006; Decressi n et al.l 
2007b,a). Gratton et al. (2001) and Carretta et al. (2005) uncov- 
ered rather large variations in both O and Na, anti-correlated 
with each other, in SGB and dwarf stars in N GC 6397 (see also 
iNorris & Da Costalll995l;ICastilho et al.ll2000i for bright giants). 
This cluster has thus suffered from internal pollution in its in- 
fancy. 

In the present study we analyse Li and Na abundances for a 
large sample of stars in NGC 6397, in an attempt to disentan- 
gle the primordial value of Li, with effects of atomic diffusion 
and intrinsic stellar depletion on the one hand, and early cluster 
pollution on the other hand. § 2 describes the observations and 
data reduction, § 3 describes the determination of stellar param- 
eters and the abundance analysis. In § 4 we present the Li and 
Na abundances found, and in § 5 we discuss signatures of Li 
depletion in the cluster § 6 summarises our conclusions. 

2. Observations 

The observations include spectroscopic data, described in § 2.1, 
and photometric data, described in § 2.2. All targets, coordinates, 
and photometry are listed in Table 2. The locations of the targets 
in the observed colour-magnitude diagram (V - (v - y)) of the 
cluster are shown in the right-hand panel of Fig.|5] 

2. 1 . High and medium-high resolution spectroscopy 

All our targets have been selected from the Stromgren uvby 
photometric survey caiTied out by F. Grundahl, matching the 
(b - y) and ci ranges spann ed by the field stars analysed by 
ICharbonnel & PrimasI (|2005|) . In total, 349 stars were observed 
across the colour-magnitude diagram of NGC 6397, from just 
below the cluster TO point (V ~ 17) to the end of the red gi- 
ant branch (RGB) (V ~ 1 1.5). Each evolutionary phase is well 
sampled with around 180 stars at the TO, 80 on the SGB, and 90 
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Stars on the RGB. Here we presen t data from two obser ving runs 
with FLAMES on the VLT-UT2 ('Pasquin i et ani200 2'), one that 
collected the Li data for the whole sample in June 2006 and a 
second run that completed the data set with Na abundance indi- 
cators for a subset of the targets in Aug. 2008. The dates, set-ups, 
exposure times, and atmospheric conditions are summarised in 
Table n 

Both sets of observations have used the GIRAFFE -hUVES 
combined mode of FLAMES, with Medusa fibres allocated to 
the faintest targets and UVES fibres allocated to the brightest 
end of the RGB. During the first observing campaign (2006), 
we selected the high-resolution GIRAFFE set-ups H679.7 (also 
known as HR15) and H627.3 (HR13), which yield respectively 
7? = 19 000 and 22 300, and spectral ranges of 660.7 - 696.5 nm 
and 612.0 - 640.5 nm. Both settings were observed in com- 
bination with the UVES Red 580 setting (R = 47 000, spec- 
tral range: 480 - 680 nm, with a small gap in the middle). 
The first GIRAFFE set-up covers the ^Li resonance line at 
670.7 nm, whereas the second setting covers the Nal doublet at 
615.4/616.0nm,. 

Unfortunately, the detection of the Nal doublet at 
615. 4/616. Onm proved to be realistic only for the brightest tar- 
gets of our sample, therefore preventing conclusions on the ex- 
istence of a Li-Na anti-correlation. To overcome this problem, 
we proposed to observe a stronger Nal doublet (at 818.3nm 
and 819.4nm,) and were granted 3 hours through the Director 
General Discretionary Time channel in Aug. 2008. Because of 
the limited time request we carefully selected a sub-set of our 
initial sample, mainly including the brightest TO and SGB stars 
(in total, 1 17 stars), with both high and low Li. For these obser- 
vations, we used the high-resolution GIRAFFE set-up H805.3A 
(also known as HR19A, characterised hy R - 14000 and 
spectral range 774.5-833.5 nm,) in combination with the UVES 
Red 860 set-up (R = 47 000, spectral range 760 - 1000 nm). 

Standard data reduction of the 2006 GIRAFFE observa- 
tions was performed with the Geneva Base Line Data Reduction 
Software (girBLDRS), version 1.13.1. To correct for the highly 
elevated dark current in the upper corner of the old CCD a 
carefully scaled and smoothed 2D dark- frame was subtracted 
from each science frame. The 2008 data were obtained after the 
GIRAFFE CCD upgrade that took place in Apr /May 2008 and 
contain negligible dark current. These data were reduced with 
the ESQ GIRAFFE pipeline, version 2. 5. 3. UVES data reduc- 
tions were performed with the FLAMES UVES pipeline, version 
2.9.7. 

After initial reduction followed sky-subtraction and radial- 
velocity correction. Through radial-velocity measurements we 
identified six non-members. Furthermore, we disregarded one 
star that appeared to be too metal-rich to be a cluster member, 
three stars that clearly fall off the cluster sequence, and yet three 
more stars for which the observed spectra are of too poor quality. 
All stars that were rejected are listed with comments in Table 2. 
The remaining sample consists of 32 stars with spectra obtained 
with the UVES Red 580 set-up (we hereafter refer to this set as 
the 'UVES sample') and 284 stars with spectra obtained with the 
GIRAFFE H679.7 and H627.3 set-ups. The UVES sample and 
the GIRAFFE sample have 1 1 stars in common. A sub-sample 
of eight stars were observed in the UVES Red 860 set-up and 
1 17 with the GIRAFFE H805.3A set-up. The UVES sample and 
GIRAFFE sample have two stars in common. 

The mean barycentric radial velocity obtained for all cluster 
members is 18.5 + 0.2 kms"', in good agreeme n t with Lind et al. 
(I2008h (18.1 ± 0.3 kms ') and iMilone et alJ (l2006h (18.36 + 
0.09(+0.10)kms-'). 



Table 1. FLAMES observations. 



Date 


UVES 


GIRAFFE 


Exposure 


Average 




setting 


setting 


time [s] 


seeing ["] 


2006 Jun 19 


Red 580 


H627.3 


2 X 2775 


1.08 


2006 Jun 19 


Red 580 


H679.7 


1 X 2775 


1.08 


2006 Jun 20 


Red 580 


H679.7 


2 X 2775 


1.08 


2006 Jun 20 


Red 580 


H627.3 


2 X 2775 


1.01 


2006 Jun 21 


Red 580 


H627.3 


2 X 2775 


0.75 


/OUo Jun zi 


D.^^ con 
Red 5oU 


hLo/9./ 


i X 1/ Ij 


n no 
0. /s 


2006 Jun 26 


Red 580 


H627.3 


3 X 2775 


1.16 


2006 Jun 26 


Red 580 


H679.7 


2 X 2775 


0.76 


2006 Jun 27 


Red 580 


H679.7 


1 X 2775 


0.56 


2006 Jun 27 


Red 580 


H627.3 


3 X 2775 


0.46 


2006 Jun 28 


Red 580 


H679.7 


3 X 2775 


0.69 


2008 Aug 3 


Red 860 


H805.3A 


2 X 2775 


0.94 


2008 Aug 28 


Red 860 


H805.3A 


1 X 2775 


1.17 



2.2. Stromgren photometry 

The Stromgren uvby photometry employed for this investigation 
was obtained at the Danish 1.54m telescope on La Silla, Chile, 
and is identical to the data set used by Korn et al. (2007) and Lind 
et al. (2008). We reiterate the main points here. All uvby obser- 
vations were collected during a two week run in May 1997, with 
a large number of observations of standard stars observed on the 
photometric nights. The field covered is roughly 9' in diameter 
and slightly west of the cluster centre. For the photometric reduc- 
tions we used the same programs and procedures as described in 
Grundahl et al. (1998, 1999, 2000, 2002a). The photometric zero 
points in vby have errors of approximately 0^01. 

The observed (b - y), (v - y), and C[ colours were corrected 
for the reddening of NGC 6397, using the relations = 
0.74 X Eg V, E ,.^, = 1.70 X £■/,_,., and = 0.20 x Et-y 
(ICrawf ordI 19751). We adopted a value o f 0.179 for Eb-v, follow- 
ing Anth onv-Twarog & Twarod (l2000l) . and 12.57 for the clus- 
ter distance modulus. The observed colours of our targets show 
a spread around the cluster sequence, of the order of 2.2% for 
(v - y) and larger, 3.7%, for (b - y). To reduce the uncertainty in 
the photometric colours we constru cted colour - magn itude fidu- 
cials for the cluster, as described in iKorn et all (l2007l) and Lind 
et al. (2008). The observed colour is interpolated at constant V- 
magnitude onto the fiducial sequence. The main advantage of 
using this method is that relative effective temperatures are more 
precisely determined, which in turn lowers the spread in chemi- 
cal abundances. 



3. Analysis 

For the abundance analysis we used a grid of ID, LTE, 
plane-parallel aii d spherical, MARCS model atmospheres 
(iGustafsson et al.l 12008). LTE abundances of Li and Na were 
derived from spectrum synthesis, using the Uppsala code bsyn. 
Non-LTE abundance corrections were thereafter applied. LTE 
abundances of Fe and Ca were derived from equivalent width 
measurements, using the corresponding Uppsala code eqwidth. 

3.1. Effective temperatures and surface gravities 

The cluster Mvfey-photometry was used to calculate effective 
temperatures for the whole sample. We implemented a num- 
ber of different calibrations based on the colour indices (v - y) 
and (b - y) and ci . Table |2] lists effe ctive tempera t ures c alcu- 
lated using the relations published in lAlonso et alj (Il996h and 
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Table 2. Photometry and effective temperatures. The full table can be retrieved from CDS/Vizier. 



ID 


RA (J2000) 


DEC (J2000) 


IV 


V 


(fc-.v)o 


(v - v)o 


Co 


Teff O09 


Teff RM05 


Teff A9699 


Teff Ha 








[kms-'] 










[K] 


[K] 


[K] 


[K] 


12138 


17 40 17.23 


-53 41 57.37 


13.5 


11.493 


0.604 


1.379 


0.452 


4727 


4660 


4650 


4540 


16405 


17 40 28.64 


-53 38 32.42 


13.5 


11.771 


0.582 


1.319 


0.429 


4790 


4719 


4709 


4564 


14565 


17 40 24.07 


-53 42 52.74 


17.6 


11.861 


0.575 


1.300 


0.422 


4811 


4738 


4728 


4569 


5644 


17 39 52.51 


-53 39 11.24 


16.6 


11.924 


0.570 


1.286 


0.417 


4826 


4751 


4742 


4576 



6500 
6000 
5500 
5000 
4500 



\ 

N 

\ 




RM05 

A9699 

O09 














■■.vO 









































0.8 



1.0 1.2 



1.4 0.3 



0.4 0.5 

(b-y)o 



0.6 



Fig. 1. The {b - y) - T^ff and (v - y) - T^s relations considered 
in the analysis. The vertical lines mark the area that is linearly 
interpolated (see text). The O09 scale forpre-TO stars fall above 
the one for post-TO stars and vice versa for A9699. 



^> 2 - 




0.3 



0.4 0.5 

(h-y)o 



0.6 



6500 6000 5500 5000 4500 
r.,, [K] 



Fig. 2. Left: The {b — y)o - My fiducial used to derive stellar 
parameters. Right: The three (/7-y)-based stellar-parameter sets. 



lAlonso et al.l (|1999|) ( hereafter these two papers are referred to 
as A9699), [R amirez & Melended (l2005h (hereafter RM05), and 
lOnehag et ail (.2009.) (hereafter O09). The first three are cali- 
brated on the infra-red flux method (IRFM), whereas the last is 
based on theoretical col ours from MARCS model atmospheres 
jGustafsson et al.l l2008l) . For all stars we assumed [Fe/H] - 
-2.0. The IRFM colour- Jeff relations constructed for main se- 
quence (MS) stars were calibrated on stars with \ogg ^ 3.8. We 
accordingly set a limit in absolute visual magnitude. My > 3.3, 
for which the MS calibrations are trusted. Analogously, the gi- 
ant calibrations were calibrated on stars with \ogg S 3.5 and we 
trust them for My < 2.8. This leaves a gap on the SGB where 
no calibration is suitable. In this range we interpolated linearly 
between each dwarf- and giant-calibration pair (see vertical lines 
in Fig. [Til. To each of the obtained eff'ective-temperature scales 
we calculated surface gravities from the relation between surface 



gravity, effective temperature, stellar mass, and luminosity. The 
details of the calculations are given in Lind et al. (2008). 

O09 have calculated synthetic (b - y)-colours for a grid of 
effective temperatures, surface gravities, and metallicities. We 
interpolated in the grid constructed for [Fe/H] = -2.0, which 
spans the parameter-space of the sample (see Fig. 4 in O09). The 
interpolation was done iteratively until the parameters found for 
each star were fully consistent, i.e. the Teff - \ogg pair found 
based on the observed colour are also constrained by relation to 
mass and luminosity. 

For the 32 RGB stars observed with UVES, we could also 
derive Ha-based effective temperature s, using bsyn. Stark broad- 
ening is based on the tabulations by Stehle & HutcheonI (Il999|) 
and self-broadening of hydrogen follow Barklem et al. (200i3)- 
The fitting method is automated and based on metal l ine-free re- 
gions extending up to ±5 nm from the line centre (see lLind et al.l 
2008), avoiding the line core. The wings of the Ha-line become 
very narrow high up on the RGB, but some sensitivity to effec- 
tive temperature remains also for these stars. With a high S /N- 
ratio it is possible to constrain Teff also for the coolest stars in 
our sample. 

Figure |2] shows a comparison between all {b - y) - Teif and 
(v - y) - Teff relations. The O09 (b - 3')-calibration matches the 
slope of the corresponding IRFM-based calibrations well, but is 
offset to higher effective temperatures by approximately 100 K. 
However, a recent updated calibration based on the IRFM finds 
effective temperatures that are higher than the existing IRFM 
calibrations, by similar order of magnitude, at these metallicities 
(Casagrande et al. 2009, in preparation). Interestingly, the O09 
scale predicts slightly higher effective temperatures for dwarfs 
compared to SGB stars of the same colour, whereas the IRFM 
of A9699 has the opposite behaviour (through the sensitivity to 
the ci -index). The interpolation between dwarf and giant cali- 
brations appears more suitable for the {b - y)-relations. The two 
(v-y) MS scales are too flat to match the slope of the correspond- 
ing giant relations well. For this reason we regard the (b - y)- 
relations as more trustworthy and place greater weight on those. 
The Hff-based values are generally cooler than all photometric 
calibrations, but the relative agreement is satisfactory. 

The right panel of Fig. |2] shows all (b - 3')-based effective- 
temperature scales with corresponding surface gravity values. 
Apparently, the two IRFM scales agree well for all but the hottest 
stars, where the A9699 scale predicts cooler effective tempera- 
tures. Again it is seen that the theoretical MARCS scale is offset 
to higher effective temperatures. The left panel of Fig |2] shows 
the cluster fiducial for (b - y)o - My- 

3.2. Metallicity 

To constrain the cluster metallicity and microturbulence val- 
ues, we calculated iron abundances from Fel and Fell lines 
for the UVES sample. The ionisation and excitation equilib- 
rium of Fe also provide an additional check of the derived 
stellar parameters. We selected 48 unblended Fel lines that 
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Table 3. Iron lines measured in the UVES spectrum of the RGB 
star #17691. The full table can be retrieved from CDS/Vizier. 



10 



Ion 


A [nm] 




Ref. 


W, [pm] 


A(Fe) 


Fel 


491.9 


-0.34 


1 


9.43 + 0.12 


5.38+0.03 


Fel 


498.6 


-1.33 


1 


4.76 + 0.05 


5.38+0.01 


Fel 


499.4 


-2.97 


1 


8.24 ±0.10 


5.55+0.03 


Fel 


501.5 


-0.30 


1 


4.00 ± 0.07 


5.41+0.01 


Fel 


502.8 


-1.12 


1 


1.97 + 0.05 


5.36+0.01 


Fel 


504.4 


-2.02 


1 


1.82 + 0.08 


5.39+0.02 
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-5.43 


1 


2.22 ± 0.09 


5.64+0.02 


Fel 


508.3 


-2.84 


1 


8.52 ± 0.09 


5.53+0.02 


Fel 
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-3.25 


1 


7.19 + 0.05 


5.56+0.01 


Fel 


514.2 


-2.24 


1 


3.35 + 0.06 


5.49+0.01 


Fel 


519.9 


-2.09 


1 


5.19 + 0.11 


5.47+0.02 


Fel 
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-1.16 


1 


3.95 + 0.08 


5.42+0.02 


Fel 


522.6 


-4.76 


1 


4.71+0.09 


5.60+0.02 


Fel 


523.3 


-0.06 


1 


10.25 + 0.10 


5.30+0.02 


Fel 


524.2 


-0.97 


1 


2.54 + 0.05 
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Fel 
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-4.97 


1 


4.03 + 0.09 
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Fel 


530.7 


-2.91 


1 
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Fel 


534.0 
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1 


5.98 + 0.09 
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Fel 


536.5 
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1 


3.48 + 0.05 
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Fel 


536.5 


-1.02 


1 


1.95 + 0.09 
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Fel 
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0.64 


1 


5.55 + 0.04 


5.16+0.01 


Fel 
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-1.98 


1 


12.90 + 0.09 
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Fel 


540.6 
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Fig. 3. Left: Histogram of the slope of Fe abundance with ex- 
citation potential of neutral and singly ionised Fe lines. Only 
the UVES sample is shown. Solid lines represent the O09 effec- 
tive temperature scale and dotted lines the A9699 scale. Right: 
Histogram of Fe abundance derived from Fe II lines minus the 
Fe abundance derived from Fe I lines. 



(l)O'Brianet 
(3) Raassen & 



al. (1991) 
Uvlings (1998) 



(2) Biemont et al. (1991) 
(4) Fuhretal. (1988) 



have oscillator strengths determined by lO' Brian et al.l (Il991h 
and 14 Fe II lines th at have oscillator strengths dete rmined by 
iBie mont et al] ( Il99ll) . iRaassen & Uvlingsl d 19981) . or lFuhret all 
(1988^ Equivalent widths were measured by fitting Gaussian 
line profiles to the observed, fixing the FWHM of lines weaker 
than 5 pm to the average value for each individual star In the 
GIRAFFE spectra only very few Fe I and Fe II lines are avail- 
able and it is not possible to obtain reliable measurements for the 
fainter half of the sample. Table[3]lists all Fe I and Fe II lines with 
references to the adopted log gf values, as well as the equivalent 
widths and abundances derived for one RGB star The full table, 
including all UVES targets, can be retrieved from the CDS. The 
table lists also the average abundance derived from Fe I and Fe II 
lines, respectively, with the error bar representing the standard 
deviation of the mean abundance. 

Microturbulence (xt) values were determined by requiring 
that Fe abundances show no trend with reduced equivalent width 
(log(W^//l)) of Fe I lines. We empirically constructed a linear pa- 
rameterisation of microturbulence with surface gravity, using the 
UVES targets only, and adopted it for the whole sample. This 
involved extrapolation to the SGB stars and dwarfs, since our 
observations do not cover enough strong lines that are sensitive 
to microturbulence in these stars. The impact on the inferred Li 
abundances is however negligible. The largest equivalent width 
measured for the Li line is 6.18pm, for the SGB star #13160. 
Varying the microturbulence between 1 km s"' and 2 km s"' low- 
ers the Li abundance inferred for this star by only 0.01 dex. 

The excitation equilibrium of Fe is commonly used to con- 
strain effective temperatures. However, 3D modelling of stellar 
atmospheres show that especially low excitation lines may give 
stron gly overestimat ed abundanc es in ID mo dels of metal-poor 
stars lAspIund et"an[T99 9; Colle t et al.ll2007l) . In addition, non- 
LTE effects are likely to ha ve diff erent impact on lines of differ- 
ent excitation potential (e . g. I Asplund 2005 ) . The extent to which 
this could bias the excitation equilibrium is presently not clear 
The left-hand histogram in Fig. [3] shows the slope of Fe abun- 
dances with excitation potential of Fel lines, derived for the 
A9699 and O09 scales. Apparently, both scales are too hot to 
satisfy the equilibrium. The necessary corrections are of order 
(-100)- (-200) K. 

The right-hand histogram in Fig [3] shows the Fe abun- 
dance derived for singly ionised Fe lines minus the Fe abun- 
dance derived from neutral Fe lines. The A9699 scale shows 
an offset from zero of approximately -0.14 dex, whereas the 



6 



K.Lind et al.: Intrinsic Li depletion and Li-Na anti-correlation in NGC 6397 



O09 scale has a smaller offset, -0.05 dex. The mean Fe abun- 
dance, as determined from Fe II lines, is 5.41 ± 0.01, using both 
scales, with a Icr dispersion of 0.04 dex. Adopting a solar abun- 
dance A(Fe)Q - 7.51 (Asplund et al. 2009 in prep.), we ob- 
tain [Fe/H]= -2.10. This is in good agreement with most re- 
cent estimates made for RGB stars in this cluster; Korn et al 
2007) find -2.12 + 0.03 for their RGB sample, Gratton et al." 
200 ll) [Fe/H]= -2.03 + 0.02 + 0.04, and Castilho et al. (2000) 
[Fe/H]= -2.0 + 0.05. Lind et al. (2008) report a lower value, 
[Fe/H]si -2.3, for stars at the base of RGB. The last study im- 
plements a cooler eff'ective temperature-scale and base the abun- 
dances mainly on Fe I lines, which explains the offset. 

The investigation is not indicating that one of the stellar- 
parameter sets is clearly preferable to the other, but we used the 
O09 scale to derive abundances of all other elements, because 
this is the most homogeneous. The changes in Li abundances 
when implementing the IRFM stellar-parameter scales are dis- 
cussed in § 5.2. 



3.3. Lithium 

The ^Li resonance line at 670.7 nm has two fine-structure com- 
ponents, separated only by 0.015 nm. At the resolution of 
GIRAFFE (R - 19300) they are indistinguishable, but a small 
line-asymmetry can be seen in the higher-resolution UVES spec- 
tra (R = 47000). Accurate measurements of the line require 
spectrum synthesis (see example fits of the line in FigHJi. We 
adopted logg/ = -0.00 9 and -0.309 for the two c omponents 
respectively, according to lLindgard & NielsonI (Il977l) . 

To derive Li LTE abundances for the stars we first nor- 
malised each spectrum by iterative fitting of a first order poly- 
nomial, taking into consideration a wavelength region extend- 
ing ±0.2 nm from the line centre. Synthetic line profiles were 
convolved with a Gaussian profile to account for the instrumen- 
tal broadening and broadening due to rotation and macroturbu- 
lence. The width of the convolving profile was fixed for each star 
by measurements of single lines of other species. The Li abun- 
dance was then varied to find the best fitting synthetic profile for 
each star through ;^f^-minimisation. In total, Li abundances or up- 
per abundance hmi ts co uld be obtained for 30 5 stars. Following 
iNorris et"aL I (l200lh and lBarklem etal.l (l2005h . errors due to ob- 
servational uncertainties were estimated with the expression: 
CTjy^ = A ■\f{n)/(R X S /N), where n denotes the number of pix- 
els spanning the full-width at half maximum for the line, R the 
spectral resolution, and S /N the signal-to-noise ratio per pixel. 

Non-LTE abu ndan c e cor rections for Li were calculated ac- 
cording to Lind et alj (1200 9), a study that is based on cal- 
culations on the same MARCS model atmospheres as here 
dGusteftson et al.' '2008^ In contrast to earlier non-LTE anal- 
yses (e.g. Carlsson et al] Il994t iTakeda & Kawano moto 2005) 
the study makes use of rigorous quantum mechanical calcula- 
tions (as opposed to the debated but common classical recipes) 
of cross-sections for collisions with neutral hydrogen, which 
have influence over the statistical equilibrium of Li. Lind et 
al. (2009) show that including the charge transfer reactions be- 
tween Li and hydrogen in the non-LTE calculations gives abun- 
dance corrections that are lower by almost -0. 1 dex for dwarfs at 
[Fe/H]= -2.0, as compared to neglecting them. For our sample, 
the MS, TO, and early SGB stars all have non-LTE corrections 
that are very similar, around -0.06 dex. With decreasing effec- 
tive temperature and surface gravity the corrections change sign 
and reach a maximum of +0.13 dex for the coolest RGB star 



3.4. Sodium 

Sodium abundances or upper abundance limits were derived 
from the 2008 observations for a subset of 117 stars, mainly 
TO stars and early SGB stars. The abundances were based on 
the Nal 818.3nm and 819.4nm doublet (when the former was 
too weak to be detected, only the latter was used). This doublet 
is second in strength after the resonance doublet at 588.5 nm 
and 589.0 nm, which is strongly affected by interstellar extinc- 
tion for the cluster stars and is therefore not useful as abun- 
dance indicator As mentioned above, the 568.2/568.8 nm and 
615.4/616. Onm doublets are too weak to be suitable for abun- 
dance analysis in the TO region at this metallicity. Despite its 
unfortunate location in the middle of a strong telluric band, the 
818. 3/819. 4nm Nal doublet does not suffer strongly from atmo- 
spheric blends that would bias the abundance analysis. 

The line-fitting method is the same as described for 
Li and oscillator strengths by iKurucz & PevtremannI (Il975h 
were adopted. We apply non-LTE abundanc e corrections to 
the 81 8. 3/819. 4nm Nal lines according to Mashonkina et all 
(|2000).The size of the corrections at this metallicity ranges from 
-0.14 dex for the hottest, highest surface gravity stars, to approx- 
imately -0.31 dex for the coolest giants. 

3.5. Calcium 

For the complete sample of 305 stars, only the equivalent widths 
of Ca lines 612. 2nm and 616. 2nm could be reliably measured. 
We based the Ca abundance on these two lines for 305 stars, 
adopting oscillator strengths from lSmith & O'Neilll d 19751) 

4. Results 

4.1. Li abundances 

Table |4] lists the adopted stellar parameters, equivalent widths, 
and abundances. Fig. |5] shows the Li abundances inferred side- 
by-side with a colour-magnitude diagram where the spectro- 
scopic targets are marked. In agreement with the expectations 
(see § 1), there is a rather well defined abundance plateau 
among the least evolved stars, followed by a drastic drop in 
Li abundance in the middle of the SGB (My ~ 3.3). This pat- 
tern is caused by the so-called first dredge-up, i.e., the dilution 
of the external convective stellar layers with deeper hydrogen- 
processed material. A second steep drop occurs at My ~ 0.0, 
which corresponds to the luminosity of the RGB bump. A sim- 
ilar decrease in Li abundance have been identi fied in stars lo- 
cated around the RGB bump in NGC 6752 (Grundahl et dj 
2002a) as weU as in stars in the halo field (Gratton et al. 2000|. 
iCharbonnel & ZahnI (120071) describe how this second Li abun- 
dance drop can be explained by a mixing process called thermo- 
haline convection, which becomes efficient when the hydrogen- 
burning shell crosses the chemical discontinuity left behind by 
the first dredge-up, and which rapidly transports surface Li down 
to internal hotter regions where this fragile element burns. 

The average Li abundance on the plateau, i.e. only including 
stars that have not undergone dilution due to the first dredge- 
up (My > 3.3), is 2.25 ± 0.01, where the error is the standard 
deviation of the mean. The Icr dispersion is 0.09 dex. 

4.2. Lithium data from the ESO archive 

To establish a view of the Li content in NGC 6397 that is as 
complete as possible, we have searched the ESO archive for ob- 
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Fig. 4. Example fits of the Lil 670. 7nm Une and the Nal 819.4 nm line. In each panel is indicated the identification number of the 
star, stellar parameters, and LTE abundance inferred from the given Une, according to the syntax Tes/logg/A{x)/e{x), where x is 
the corresponding element and e(x) the error in abundance. Top: A Li and Na-normal TO star. Second from the top: A Li-poor and 
Na-rich TO star. Second from the bottom: GIRAFFE giant spectrum. Bottom: UVES spectrum of same giant. 

Table 4. Adopted stellar parameters, equivalent widths, and Li, Na, and Ca abundances. The full table can be retrieved from 
CDSA^izier. 



ID 


Mv 


[K] 


log? 


[km/s] 


log L/Lq 


670.1 
[pm] 


A(Li) 


818.3 
[pm] 


819.5 
[pm] 


A(Na) 


W,, 612.2 
[pm] 


W,, 616.2 
[pm] 


A(Ca) 


12138 


-1.077 


4727 


1.52 


1.61 


2.49 


< 0.44 


< 0.08 


8.85 
±0.31 


11.49 
±0.17 


4.29 
±0.03 


10.91 
±0.08 


12.19 
±0.06 


4.70 
±0.01 


16405 


-0.799 


4790 


1.66 


1.58 


2.37 


< 0.42 


< 0.17 








10.41 
±0.06 


11.64 
±0.02 


4.64 
±0.01 


14565 


-0.709 


4811 


171 


1.58 


2.33 


< 0.43 


< 0.19 








10.30 
±0.02 


11.39 
±0.07 


4.62 
±0.01 


5644 


-0.646 


4826 


1.74 


1.57 


2.30 


< 0.55 


< 0.25 








10.25 
±0.08 


11.40 
±0.09 


4.62 
±0.01 


17163 


-0.439 


4873 


1.84 


1.56 


2.21 


0.43 
±0.13 


0.19 
±0.14 








9.97 
±0.07 


11.07 
±1.2 


4.62 
±0.01 
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Fig. 5. Le/f; The abscissa shows non-LTE Li abundances infeiTed for our sample (the archival data are not included) and the ordinate 
absolute visual magnitude (My - Y - 12.57). Arrows mark Li upper limits. Right: The spectroscopic targets marked with black 
filled circles in the colour-magnitude diagram of NGC 6397. The two horizontal dashed lines mark the location of rapid decrease in 
Li abundance caused by stellar evolution (see text). 



servations of more MS, TO, and SGB stars. Especially, to put 
constraints on the size and significance of a small increase in 
A(Li) located just before the onset of the first dredge-up (see 
§ 5.1 and Fig. |6]l we would benefi t from having more stars. In 
previous analyses, iBonifacio et al.l (|2002|) analyse Li in twelve 
TO stars (of which seven stars wer e prev iously analysed by 
Irhevenin et al.]2001 a nd three stars bv lPasquini & Molaro 1996) 
and lKorn et al.l (l2007h analyse five TO stars and two SGB stars. 
We did not incorporate those stars in our analysis since it would 
only slightly have increased the number statistics. However, we 
have retrieved and analysed a set of FLAMES archival data 
from 2007 (079.D-0399(A), P. 1. Gonzalez-Hernandez). These 
include observations in the GIRAFFE high-resolution set-up 
H665.0 (HR15N, R = 17000), covering the Lil 670.7 nm Une for 
a sample of 80 dwarfs (My -4.9) and 88 SGB stars (My ~ 3.4) 
in NGC 6397. The reduction of the data was performed with 
the girBLDRS pipeline. In total, 38 MS stars and 55 SGB stars 
were identified in our uvby photometric catalogue. For the re- 
maining stars, we add a zero-point correction of 0.09 magnitudes 
to the visual magnitude scale listed in the archival data, esti- 
mated from the targets that were cross-matched between the 
catalogues. Applying the correction, the two magnitude scales 
match to typically within 0.02 magnitudes for individual targets. 
After establishing the visual magnitudes, colours are assigned to 
each target by interpolating V onto the fiducial cluster sequences, 
as was done for the main sample (see § 2.2). 

The results of the analysis of the archival data are shown to- 
gether with our own data set in Fig. |6] One more very Li-poor 
star, marked with an arrow, was identified in the archive sample. 



for which only an upper limit to Li could be derived. Especially, 
the region around My - 3.2 - 3.5 is now much better sampled, 
with the archival data confirming the presence of the small abun- 
dance increase just before the steep drop. In addition, the plateau 
is extended to fainter magnitudes. My x 5. Including this archive 
sample we have Li abundances for 454 stars in the cluster, of 
which 346 have My > 3.3. Overall the quality of the data sets 
are similar and the abundances agree very well between the sam- 
ples. In total, 14 SGB stars overlap and thus have observations 
of the Li line in the H665.0 as well as in the H679.7 setting. 
The mean difference in Li abundance between these 14 stars is 
-0.02 ±0.02 dex, i.e. there is no systematic bias between the data 
sets. 

The average Li abundance on the plateau (i.e., excluding 
stars with My < 3.3), including both our data and the archival 
data, is 2.270 + 0.005. Unfortunately Na information is not avail- 
able for the archive targets. 

4.3. The Li-Na anti-correlation 

As discussed in § 1, globular-cluster samples may suffer from 
bias due to early pollution of the intra-cluster gas by a first gener- 
ation of more massive, faster evolving stars. It is thus important 
to discriminate between first-generation stars that formed out of 
pure pristine material and second-generation stars that formed 
out of polluted material. The former should have been born with 
the cosmological Li abundance and a low Na-abundance, sim- 
ilar to that of field stars of comparable overall metallicity. The 
latter may have inherited a lower Li abundance at birth and will 
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Fig. 6. Lithium abundance against visual magnitude. The targets 
of our observations are plotted with filled circles and the archival 
data are plotted with open squares. Information about evolution- 
ary status and eff'ective temperature for various My intervals is 
also given. 



also exhibit other signatures of H-burning at high temperature, 
in particular enhanced Na abundances resulting from the NeNa 
chain in the polluter stars. 

In Fig. |7]we report the Li and Na abundances for 100 stars 
in our sample with My > 3.3, i.e., those that have not yet un- 
dergone Li dilution via the first dredge-up. Note that it is the 
first time that these quantities are determined simultaneously for 
such a large sample of TO and SGB stars in a globular cluster 
This diagram clearly emphasises the extent of the star-to-star Na 
abundance variations (by up to 1 d ex) within NGC 6397 . This is 
in agreement with the findings by ICarretta et alj (l2005l) . which 
were based on six dwarfs and nine SGB stars. From the appear- 
ance of Fig.|7] we deduce that the most Li-poor targets indeed 
show the most elevated Na levels, implying an anti-correlation 
between the abundances of Na and Li. The spectrum of the most 
Li-poor and second most Na-rich star in the sample, #6685, is 
shown in Fig. H] The abundance pattern displayed by this star 
and by at least three more stars, agrees well with the globu- 
lar cluster self-enrichment scenario and we conclude that they 
have likely been formed from polluted gas. Field stars with sim- 
ilar metallicity to NGC 6397 have typical Na ab undances in the 
range A(Na) = 3.6 - 3.9 ( Andrie vskv et al.ll200 7). It is thus rea- 
sonable to regard stars with A(Na) < 3.9 as belonging to a first 
generation of non-polluted objects. 

To find the statistical significance of the indicated Li-Na anti- 
correlation we performed linear regression between A^(Li)/A^(H) 
and A^(Na)/A^(H), taking into account measurement errors in 
both quantities as well as up per limits on Na (IDL-routine 
linmix_err , see lKellyll2007h . The analysis gave a correlation 
coefficient of -0.6 between the Na and Li abundances, which for 
a sample of this size is highly significant. The probability of two 
uncorrected variables producing such a correlation coefficient 
is less than 0.05%. However, as can be realised from Fig.|7] the 
significance of the anti-correlation is dependent on the most Na- 
enhanced stars. In fact, performing the same linear regression for 
stars with A(Na) < 4.1 no significant anti-correlation was found. 
The Li abundances among the plateau stars are thus not much 
affected by high degrees of pollution. 
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Fig. 7. Non-LTE Li and Na abundances for the sub-sample of 
100 dwarfs and SGB stars with observations in the H805.3A set- 
ting and My > 3 . 3 . The arrows mark upper limits to the Na abun- 
dance. If the effective temperature of a target has been underesti- 
mated by 100 K, the corresponding point in this diagram should 
be shifted approximately in the direction of the tilted arrow and 
in its inverse direction for overestimated effective temperatures. 

In Fig. Q it is also indicated with a tilted arrow approxi- 
mately how the location of a star would be affected by an er- 
roneous effective temperature. An error of AFeff = +100 K cor- 
responds to AA(Li) = +0.07 and AA(Na) = +0.04. EiTors in 
effective temperature therefore weakly correlate the abundances 
with each other, rather than anti-correlate them. Hence, it is clear 
that the anti-correlation itself is not an artifact from uncertainties 
in the stellar-parameter determination. Possibly even a tighter 
anti-correlation exists for all targets, but is being distorted by the 
tendency of effective-temperature errors to align the abundances. 
As the Li abundance as well as the Na abundance are based on 
lines from neutral atoms, the effective temperature is the stellar 
parameter with by far the greatest influence over the results. The 
uncertainties stemming from surface gravity, microturbulence, 
and metallicity are all negligible in the context. 

The extent of the Li and Na variations found in the presen t 
study is sim ilar to that found in NG C 6752 jPasquini et al.ll2005l) 
and 47 Tuc dBonifacio et al.ll2"007bh . and thus appears to be inde- 
pendent of the metallicity. However, the Li-Na anti-correlations 
in NGC 6752 and 47 Tuc are based respectively on nine and four 
stars only, limiting the comparison. A more detailed discussion 
of the implications for the self-enrichment of NGC 6397 will be 
performed in a separate paper where O, Na, Mg, and Al abun- 
dance determinations will also be presented for a subset of RGB 
stars in our sample. In § 5 we use the information about Na to 
constrain intrinsic stellar Li depletion processes. 

5. Discussion 

5A. Signatures of intrinsic litliium depletion 

A proposed explanation for the difference between the Li abun- 
dance found in the metal-poor halo and the BBNS prediction 
is intrinsic Li depletion due to the combined effect of gravita- 
tional settling and weak turbulence, the nature of which is still 
a matter of debate. In the case of NGC 6397, the stars on the 
Spite plateau would have to have been depleted by typically 
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Fig. 8. Group-averaged Li abundances against visual magnitude. 
The full line is drawn for the whole sample (including the 
archival data), i.e., for both first-generation Na-poor stars and 
second-generation Na-rich stars. Each bin contains 40 stars. The 
dashed and dotted lines correspond to diff'erent sub-sets of our 
own sample, selected according to their corresponding Na abun- 
dances. Each bin contains 13 stars. In the upper left corner mean 
error bars are shown for each line. Note that the range between 
Mv « 4.4 - 4.8 is not sampled by the data. 
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Fig. 9. Group-averaged Li abundances against eff'ective temper- 
ature. The lines correspond to the same sub-samples as in Fig. [8] 
Filled circles represent post-TO stars and open squares pre-TO 
stars. In the upper left corner mean error bars are shown for each 
line. 

AA(Li) = 2.72 - 2.27 - 0.45 dex, corresponding to almost a 
factor of three. To constrain the physics involved in producing 
such depletion, it is important to accurately establish trends of Li 
abundance with evolutionary phase and effective temperature. In 
the following, we carry out our analysis based on sub-samples of 
data defined according to Na content. As discussed in § 4.3, it is 
reasonable to consider that the stars born out of pristine material 
are those with A(Na) < 3.9. 

Fig. |8] shows Li abundance vs visual magnitude on the 
plateau, with abundances binned in groups of targets for clar- 
ity. In the lower panel of each plot, the Icr dispersions in the 
corresponding bins are shown. The full lines are drawn for the 
whole sample including both our data and the archival data, i.e., 
it does not discriminate between first-generation Na-poor stars 
and second-generation Na-rich objects. The single plateau-star 



for which only an upper limit could be inferred is not included 
in the binning. The dashed lines correspond to first-generation 
stars with A(Na) < 3.9 and dotted lines to second-generation 
stars with A(Na) > 3.9. Mean error bars for each line are also in- 
dicated in the figure. Fig.|9]shows the same plot, but using effec- 
tive temperature as the reference scale instead of visual magni- 
tude. To separate pre- and post-TO stars with the same effective 
temperature, the dwarfs are represented by open square symbols 
and the SGB stars with filled circles. The approximate evolution- 
ary status of the stars can be recovered from Fig. |6] Stars with 
My > 4 are dwarf stars below the TO, having effective tempera- 
tures reaching from the maximum 6400 K to 6200 K at My ~ 5. 
Stars with My < 4 are slightly more massive objects that have 
passed the TO, i.e., SGB stars with T^ff x 6400 - 5900 K. At 
My ~ 3.3 the Li dilution process through the first dredge-up sets 
in. 

As can be seen in Figs [8] and |9] trends of Li abundance 
with evolutionary phase and effective temperature are not sig- 
nificantly different between first- and second-generation objects. 
The mean value of A(Li) appears generally to be slightly higher 
for the first-generation stars (dashed lines) than for the second 
generation (dotted lines), in agreement with what is expected in 
the globular cluster self-enrichment framework, but the trends 
are fully compatible within the error bars. We conclude that the 
mean Li abundance trends of the full sample are not significantly 
biased by intra-cluster pollution in NGC 6397. On examination 
of Fig [8] we especially note the presence of an upturn in A(Li) 
just before the steep abundance drop at My ^ 3.3, which ap- 
pears rather robust against varying A(Na) in the sample. In a 
very limited magnitude-range the mean Li abundance increases 
by around 0. 1 dex up to A(Li) - 2.35. As mentioned previously, 
a difference of this size between TO and SGB stars in NGC 6397 
was identified by Korn et al. (2007). 

The behaviour of Li abundance with T^ff is illustrated in Fig. 
|9] There is an overall increase with decreasing effective temper- 
ature in the range 6400 - 6100 K, where a maximum is reached, 
followed by a decrease again until Teff a: 5900 K, where the 
first dredge-up sets in. The initial increase appears to be present 
both for pre- and post-TO stars. Speculatively, the lower value 
found for hotter dwarfs compared to cooler dwarfs, could be 
interpreted as the counterpart of the very right-shoulder of the 
so -called Li-dip observ ed in Population 1 stars (see references 
in lTalon & Ch arbonnel 1998, 2003). This feature appears in the 
present data at exactly the same effective temperature as in open 
clusters: In the Hyades, the mean Li abundance decreases by 
approximately 0.10-0.15 dex in the effective temperature range 
between 62 00 and 6400 K. As discussed by Talon & Charbonnel 
(Il998ll2004l) . a similarity between Pop I and Pop II dwarfs is the- 
oretically expected, since the depth of the stellar convective en- 
velope, and thus the nature and the efficiency of the Li depletion 
mechanism, are expected to depend only on the stellar effective 
temperatur^ If this feature is real (see § 5.2), it seems that we 
may have discovered for the first time in Population II stars the 
very beginning (in terms of effective temperature) of the Li-dip, 
which may be more shallow at this low metallicity. Note that one 
should be cautious with interpretations of Li trends for dwarfs in 
our sample. As mentioned in § 2.1, due to observational limi- 
tations the majority of our sample stars with Na determination 



- At a given effective temperature, Pop II stars along the plateau have 
lower masses than Pop I stars on the cool side of the Li dip. Regarding 
the depth of the convective envelope however, their lower metallicity 
compensates for the mass effect. See iTalon & Charbonnell ( |2004) for 
more details. 
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are SGB stars, only a handful being relatively warm dwarfs. No 
Na abundance determination has thus been carried out for the 
coolest dwarfs below the TO. However, we have seen that the 
bias in Li abundance trends among SGB stars due to pollution is 
weak, which is likely true also for dwarfs. 

In § 4 we concluded that the most Li-deficient stars in our 
sample likely belong to a second generation of stars, having 
experienced intra-cluster pollution. Among the first-generation 
stars with low Na levels, there are no stars with Li abundances 
lower than A(Li) - 2.0. The lower panels of Figs. [8] and |9] 
show that the abundance dispersion of the full sample is always 
rather low, below 0.1 dex, except for the bin at My a; 3.7 and 
Teff ~ 6350 K, which contains the most Li -poor star with an 
abundance detection (see Fig. |6]l and therefore has a higher dis- 
persion. Focusing on the first generation of Na-poor stars only, 
the typical abundance scatter is 0.09 dex, not following any ob- 
vious trends with effective temperature or visual magnitude. Not 
shown in the figures is the typical measurement error in abun- 
dance, stemming from photon noise, which ranges from 0.08 dex 
for the hottest stars to 0.05 dex for the coolest plateau stars. To 
this measurement error should be added the propagated uncer- 
tainty in stellar parameters. Assuming that the typical star-to- 
star error in effective temperature, which is the most influential 
parameter, lies in the range 50 - lOOK, a corresponding addi- 
tional spread of 0.04 - 0.07 dex in Li abundance is expected. 
Given these basic estimates, we conclude that the observations 
are compatible with zero scatter in Li abundance among the first- 
generation stars. 

5.2. Effects of stellar parameters arid rion-LTE 

We now briefly discuss the influence on the abundance trends 
presented in § 5.1 from the choice of effective-temperature scale 
and from the non-LTE treatment. Basically, both eff'ects have a 
systematic influence over the Li data. If we were to adopt the 
IRFM-based effective-temperature scale by RM05 (see § 3.1) 
the Li abundances of post-TO stars would be systematically low- 
ered by -0.06 dex, thus aggravating further the overall discrep- 
ancy to the primordial value. The A9699 scale would act in the 
same direction, but also have a minor differential impact on the 
Li abundances since this effective-temperature scale is steeper 
in the range from the TO to further up the SGB. The TO star 
abundances would therefore be lowered by a greater amount 
(-0. 1 dex) than the cooler SGB star abundances, which in turn 
would enhance the abundance upturn seen in Fig. [8] and the cor- 
responding maximum abundance seen in Fig.|9]at T^ff - 6100 K. 
Nothing would thus qualitatively change in the discussion of 
abundance trends for post-TO stars. However, implementing ei- 
ther of the two IRFM-based scales would lower the abundances 
found for cool dwarfs more than hot dwarfs. The difference be- 
tween these two groups discussed in § 5.1, which we speculated 
could be a signature of the Li-dip, would hence be erased. 

As mentioned in § 3.3, the non-LTE abundance corrections 
are similar, approximately -0.06 dex, for all plateau stars. The 
LTE abundances therefore have the same relative behaviour but 
are offset to a slightly higher values. 



5.3. Comparison to other studies 

Comparing the found Li abundances to other determinations 
made for NGC 6397, we find good agreement with Korn et 
al. (2007), who place their five TO stars with My ~ 3.6 at 
A(Li) = 2.25 and two SGB stars with My ~ 3.3 at A(Li) = 2.36. 
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Fig. 10. Top: Comparison between bin-averaged Li abundances 
(red filled circles connected with solid lines) and the predictions 
from the stellar-structure models of Richard et al. (2005). T5.80 
represents the model with lowest efliciency of turbulent trans- 
port, T6.00 intermediate efficiency, and T6.09 highest efficiency. 
The reference scale is logarithmic luminosities in units of so- 
lar luminosities. Bottom: The same plot for Ca abundances. A 
colour version of this figure is available in the online edition of 
the journal. 



The similarity between the abundances for TO and SGB stars 
is partly due to the cancellation between the effects from the 
cooler effective-temperature scale of Korn et al. and the fact 
that they do not correct for non-LTE. For RGB stars the two ef- 
fects do not cancel (non-LTE corrections are positive for RGB 
stars) and Korn et al. thus find abundances that are approxi- 
mately 0.1 dex lower than ours. The Li abundances reported by 
Bonifacio et al. (2002) for TO stars lie on average 0. 1 dex higher 
than ours. This difference can be fully traced to their effective- 
temperature scale, which is hotter by approximately 50 K than 
the one adopted in this study, and to t heir use of different non- 
LTE corrections CCarlsson et alJ[T994b . For these stars, the cor- 
rections by Ca rlsson et aL are sm aller in absolute value than the 
ones found bv lLind etal] (|2009^. The Li abundances and abun- 
dance upper limits reported by Castilho et al. ( 2000) for 1 6 RGB 
stars are in good agreement with ours, except for the Li detec- 
tions made for giants cooler than 4900 K, which are significantly 
higher than what we find. 

In the extensive literature study of field stars carried out 
bv lCharbonnel & PrimasI (l2005b . the authors concluded that the 
stars in their TO and SGB samples generally are more Li- 
rich than their dwarf sample. Especially, dwarfs that are cooler 
than 6000 K show a distinct decrease in Li abundance with de- 
creasing effective temperature. As the dwarfs in our sample are 
all hotter than 6000 K we cannot verify this result. However, 
with hindsight one may trace a tendency of increasing Li abun- 
dance with decreasing effective temperature in the range T^ff = 
6050 - 6400 K among the MS, TO, and post-TO stars observed 
by Charbonnel & Primas (see their online tables 5-9), similar to 
the one we find. 
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Fig. 11. Top: Same as in in Fig[TO]but with effective temperature 
as the reference scale. Here, only post-TO stars are included. 
Bottom: The same plot for Ca abundances. A colour version of 
this figure is available in the online edition of the journal. 



5.4. Comparison to diffusion-turbulence models 

Fig. [TOland Fig. [TTI show comparisons of our obtained Li abun- 
dances with the predictions from stellar-structure models in- 
cluding atomic diffusion fr om first principles and an ad-hoc 
recipe for turbulent mixing dRichard et al .1120051 and references 
therein). For clarity we show the same group-averaged Li abun- 
dances as in Fig.[8]and Fig.|9l including the whole sample, irre- 
spectively of Na abundance. Only post-TO stars appear in Fig. 
nn For comparison, the bottom panels of both figures show cor- 
responding trends for binned Ca abundances. Visual magnitudes 
are converted^ to stellar luminosities using the bolometric cor- 
rections of lAlonso et al.l (Il999i) . The same three models, T5.80, 
T6.00, and T6.09, with different efficiencies of turbulent trans- 
port, as displayed in Korn et al. (2007) and Lind et al. (2008) 
are shown. To ensure that the model Teff - log(L/L0)-plane is in 
agreement with our study regarding the location of the TO-point, 
we add -nSO K to the effective-temperature scale of Richard et al. 
and -1-0.05 dex to the logarithmic luminosity-scale. The absolute 
abundance scale of the models is shifted to agree with the obser- 
vations for stars evolved beyond the onset of the first dredge-up. 

The efficiency of turbulent mixing determines how much Li 
that is transported between the external convection zone to the 
region in the stellar interior where Li is destroyed, which in turn 
has influences over the appearance of the Li abundance plateau. 
In the T6.09 model, with highest efficiency, the plateau is flat, 
whereas the surface abundances predicted by the T5.80 model 
show a clear dependence on stellar luminosity and effective tem- 
perature. As seen in Fig. [TOl the Li abundances of NGC6397 
strictly limit the efficiency of turbulent mixing, although no 
model perfectly reproduces the observed abundances. In terms 
of the size of the Li variations on the plateau, the T6.00 model 
appears to be the best choice, supporting the findings by Korn 
et al. (2007) and Lind et al. (2008). The observed Ca abundance 
trend is in good agreement with both higher-efficiency models 
(T6.00 and T6.09), and especially the location of the upturn at 
log{L/LQ) ~ 0.65 is matched. However, Fig. [TTI shows clearly 



that the behaviour of the post-TO Li abundance with effective 
temperature is not reproduced by any of the models. Especially 
the location of the maximum A(Li) occurs at higher effective 
temperature than the models predict. 

The initial Li abundance predicted by the models, account- 
ing for the adjustment we make in the vertical direction, is 
A(Li) - 2.46, which is significantly closer, but still not in full 
agreement with the latest predictions from BBNS of 2.72 + 0.06 
jCvburt et alj|2008h . 



6. Conclusions 

We have presented a comprehensive study of the Li content in 
a large sample, 454 stars in total, of MS, TO, SGB, and RGB 
stars in the metal-poor globular cluster NGC 6397. The cluster 
dwarfs and early SGB stars form a thin Li abundance plateau on 
the same level as the Spite plateau for field stars, whereas more 
evolved stars have undergone drastic Li depletion due to phys- 
ical processes connected to low-mass stellar evolution. Two re- 
markably well-defined locations in absolute visual magnitudes, 
Mv ~ 3.3 and My ~ 0.0 respectively, are identified as the lo- 
cations of corresponding steep Li abundance drops. They corre- 
spond respectively to the occurrence of the first dredge-up on the 
SGB and to the onset of th ermohaline mixing at the RGB bump 
jCharbonnel & Zahnll2007l) . 

Using information of Li and Na abundance for a sub-sample 
of 100 dwarfs and early SGB, we have for the first time iden- 
tified a significant anti-correlation between the two elements in 
this cluster. This is interpreted as the signature of intra-cluster 
pollution from a previous generation of more massive stars. The 
spread in Na and Li abundances of stars not having undergone 
Li depletion due to dredge-up, is very large, about one order- 
of-magnitude. However, only a handful of stars show signifi- 
cantly depleted Li levels and the identified anti-correlation de- 
pends critically on these objects. For Na-enhancements up to 
0.7 dex, no corresponding Li-deficiency can be detected. The av- 
erage Li abundance is thus not much affected by internal pollu- 
tion. These abundance patterns have consequences for for the 
self-enrichment scenario in NGC 6397, which will be discussed 
in a forth-coming paper 

By dividing our sub-sample with Na determinations into two 
groups, consisting of first and second generation of stars ac- 
cording to the degree of pollution, we illustrate how especially 
trends of Li abundance with effective temperature and evolution- 
ary phase are in reasonable agreement between the two genera- 
tions. The average Li abundance is typically A(Li) - 2.25 for 
stars located below and above the turn-off, and show a slight up- 
turn of AA(Li) a;0.1 dex for stars in the middle of the subgiant 
branch. This difference agrees with with previous findings for 
the cluster as well as for field stars. We find no support for a 
significant Li dispersion among the first generation of stars. 

The identification of a minor Li deficiency of the hottest 
dwarfs in our sample compared to slightly cooler dwarfs leads 
us to suggest that the very right-hand wing of the Li-dip may be 
visible also among Pop 11 stars. However, the signature is erased 
if a different effective temperature scale is adopted. 

The detailed picture we have formed of the Li abundance 
trends in this globular cluster can be used to constrain the physics 
involved in depleting Li from the photospheres of low-mass 
metal-poor stars, as well as the extent of the depletion. To il- 
lustrate this, we compare our observational findings to predic- 
tions from stellar-structure models including atomic diffusion, 
with additional turbulence below the convection zone. We find 
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that some turbulence, in a very limited efficiency-range, is in- 
deed required to explain observations. However, these models 
fail to reproduce the behaviour of Li abundance with elfective 

temperature along the plateau, suggesting that a detailed under- 
standing of the physics responsible for depletion is still lacking. 
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